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Recently, metamaterials with near-zero permittivity (epsilon-near-zero, ENZ) emerge into 
the focus of the extensive explorations in both theory and engineering due to their anoma- 
lous electromagnetic features at microwave and optical frequencies. ENZ metamaterials 
have been widely used in lots of exciting applications, such as directive emission 1 , electro- 
magnetic energy squeezing and tunneling 2-5 , electromagnetic wave phase front shaping 6 ' 7 , 
electromagnetic transparency 8 ' 9 , and invisible cloaking 10 ' 11 . In the previous work, it was 
shown theoretically that in ENZ media loss can induce transparency, omni-directional col- 
limation, and transmission enhancement due to the interplay between ENZ and material 
loss 12 . Anisotropic loss can even improve the propagation of oblique incident beams 12 . In 
this Letter, we verify such abnormal transmission enhancement in anisotropic ENZ metama- 
terials with full-wave numerical simulations together with the iso-frequency contour (IFC) 
analysis, and demonstrate the possibility in experimental realization of this phenomenon by 
using silver-germanium multilayered structures. 

The propagation of a plane electromagnetic wave incident upon the flat interface (along 
the x-direction) between air and an anisotropic ENZ metamaterial is considered in two- 
dimensional space. The permittivity of air is set to be unity, and the principle components 
of the permittivity tensor for the anisotropic ENZ metamaterial are denoted as e x and s y: 
respectively. The permittivity e x is also set to be unity, i.e., e x — 1, while the permittivity 
s y reads e y = Re(£ y ) + ilm(£ y ) with a near-zero real part Re(£ y ) and a non-zero imaginary 
part Im(£ y ) (the material loss). Regarding a plane electromagnetic wave in TM polarization 
(with non-zero field components E X) E y) and H z ), the dispersion relation of the anisotropic 
ENZ metamaterial reads 

k 2 k 2 

f + f = ki (i) 

£y ^x 

in which fc is the wave vector in free space. The wave vector component k x is real and 
determined by the incident wave, while k y is complex and it is denoted as k y = Re(k y ) + 
ilm(/c y ), with the imaginary part lm(k y ) representing the propagation loss. The iso-frequency 
contours (IFCs) of the anisotropic ENZ metamaterial can be plotted based on the dispersion 
relation in Eq. (1). Figure 1 displays three different IFCs at the frequency of 193.4 THz 
(1.55 /xm) corresponding to different values of permittivity e y) which has the same real part 
of Yle{e y ) = 0.001, but different material losses: the low loss lm(s y ) = 0.002 (Fig. 1(a)), the 
moderate loss Im(£ y ) = 2 (Fig. 1(b)), and the high loss lm(s y ) = 20 (Fig. 1(c)), respectively. 
The IFC of k x ~ Re(k y ) is plotted as red solid curve, while the IFC of k x ~ lm(k y ) is 
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plotted as red dashed curve. The circular IFC of air with radius ko is also plotted as 
blue solid curve for comparison. Additionally, the incoming electromagnetic wave in air is 
represented by the wave vector k- in and the Poynting vector £ in , while the electromagnetic 
wave inside the anisotropic ENZ metamaterial is represented by the wave vector fc ENZ and 
the Poynting vector S wz . As a particular example, k- in and S'in are marked as the incoming 
electromagnetic wave in air has an incident angle of 45°. 

It is clear that for the low loss Im(£ y ) = 0.002 in Fig. 1(a), the IFC of k x ~ Re(k y ) forms 
a hyperbolic-like shape, with a very narrow opening toward the /^-direction. Therefore, the 
direction of the transmitted power inside the ENZ varies from normal to nearly parallel to 
the interface. Moreover, the high propagating loss, represented by the IFC of k x ~ lm(k y ) 
prevents the electromagnetic wave from transmitting inside the anisotropic ENZ metamate- 
rial even at small incident angles. Finite element method (FEM) full- wave simulations are 
used to illustrate the wave propagation behavior, where a TM polarized Gaussian beam is 
applied to mimic the experimental condition instead of the ideal plane electromagnetic wave. 
The Gaussian beam is represented by the distribution of the H z component. Four different 
angles of incidence are calculated, including 15°, 30°, 30°, and 45°. The simulation results 
clearly indicate that for all these angle of incidence, the incoming Gaussian beam cannot 
couple into and transmit in the anisotropic ENZ metamaterial, but reflect back at the in- 
terface, which is consistent with the above IFC analysis. Thus, a low loss of the anisotropic 
ENZ metamaterial along the ^/-direction will prevent the transmission due to impedance 
mismatch. As the material loss is increased to the moderate loss Im(£ y ) = 2 in Fig. 1(b), 
the opening of the IFC of k x ~ Re(k y ) toward the /^-direction becomes wider, leading to a 
flatter curvature. Besides, the corresponding propagating loss for the oblique incidence also 
decreases due to the increased material loss 12 . Additionally, because of the flatter curvature 
of the IFC of k x ~ Re(k y ), the incoming electromagnetic wave with small angle of incidence, 
e.g., less than 45°, will transmit inside the anisotropic ENZ metamaterial with the direction 
nearly normal to the interface, since the Poynting vector of the electromagnetic wave should 
be perpendicular to the IFC of k x ~ Re(k y ). In the simulations, the Gaussian beam can 
transmit into the anisotropic ENZ metamaterial with a wide range of the angles of incidence 
from 15° to 60°. However, it is worth noting that the propagation distance of the Gaussian 
beam inside the anisotropic ENZ metamaterial is limited by its propagation loss, especially 
for the Gaussian beam with large angle of incidence, which is indicated in the IFC. Fur- 

3 



thermore, the simulations also show that the wave propagation direction (indicated by the 
Poynting vector S ENZ ) of the Gaussian beam inside the anisotropic ENZ metamaterial is 
different from the direction of its phase velocity (indicated by the wave vector fc ENZ ), due 
to the anisotropic permittivity profile. For the high loss lm(e y ) = 20 in Fig. 1(c), the IFC 
of k x ~ Ke(ky) becomes almost flat straight lines, which allows the electromagnetic power 
inside the ENZ to propagate in the normal direction for a wide range of the incidence an- 
gle, for example, the Gaussian beam with 60° angle of incidence shown in the simulations. 
The relatively low propagation loss permits a long propagation distance of the transmit- 
ted electromagnetic wave inside the anisotropic ENZ metamaterial. As a consequence, the 
transmission is enabled and enhanced by the high loss of the anisotropic ENZ metamaterial 
along the ^/-direction. 

Based on the effective medium theory 13 , we numerically demonstrate this exotic electro- 
magnetic phenomenon with realistic metal-dielectric multilayered structures, as depicted in 
Fig. 2. The anisotropic ENZ metamaterial is constructed as a multilayer consisting of alter- 
nating thin layers of silver and germanium stacking in the x-direction. In order to reduce the 
reflection at the interface, the multilayered stack is placed on a germanium substrate. The 
thickness of each sliver-germanium pair (along the x-direction) is 100 nm, and the effective 
permittivity of the stack can be represented as 13 

£ex = (f/e A g + (1 - iO/^Ge)" 1 , 
£ey = f^Ag + (1 — /)^Gej 

where / is the filling ratio of the silver. The permittivity of the silver follows the simple 
Drude model e^g = £oo + (jJp/(cj(cj + icry)) with permittivity constant = 5.0, plasma 
frequency uj p = 1.38 x 10 16 rad/s, and damping rate of bulk silver 7 = 5.07 x 10 13 rad/s. 
The additional factor a in the simple Drude model is used to account for the increased loss 
due to surface scattering, grain boundary effects in the thin silver film and inhomogeneous 
broadening in the realistic experiment 14 ' 15 . In the simulation, three different values of factor 
a are considered, including a = 1, a = 3, and a = 9, which results in three different 
damping rates: 7 (Fig. 2(a)), 37 (Fig. 2(b)), and 97 (Fig. 2(c)). Moreover, the permittivity 
of germanium is £Qe = 19.010 + 0.087i 16 , and the ENZ frequency, at which Re(e ey ) = 0, 
is designed at 193.4 THz (1.55 /xm), thus based on Eq. (2) the filling ratio of the silver 
is / = 0.129 (Fig. 2(a)), / = 0.130 (Fig. 2(b)) and / = 0.145 (Fig. 2(c)), with respect 
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to the different damping rates. A TM polarized Gaussian beam at the ENZ frequency 
193.4 THz is considered in the simulation, and two angles of incidence, 5° and 15°, are 
calculated for each damping rate. Since lower damping rate leads to smaller material loss 
lm(e ey ) along the y-direction in the multilayered stack, it is clear that the Gaussian beam 
almost cannot couple into the silver-germanium multilayer stack and propagate inside for 
both angles of incident, although there is some optical scattering along the interface due 
to the layered structures. Furthermore, it is worth mentioning that the Gaussian beam 
splitting phenomenon shown in Fig. 2(a) inside the multilayered stack is caused by the 
optical nonlocality in ENZ metamaterials 17 ' 18 . As the damping rate increases, as shown in 
Fig. 2(b) and 2(c), more energy of the Gaussian beam can enter into the stack for both angles 
of incidence, especially for the lower angle of incidence, resulting in a clear demonstration of 
the loss enhanced transmission in realistic metamaterial structures. Moreover, the simulation 
results also indicate that the Gaussian beam propagates inside the multilayered stack normal 
to the interface between the stack and germanium substrate, which agrees well with the above 
analysis. 

To conclude, we have verified the extraordinary transmission enhancement and beam 
collimation caused by the material loss in anisotropic ENZ metamaterials through numer- 
ical simulations, and explained this phenomenon by applying the IFCs of anisotropic ENZ 
metamaterials. Furthermore, we also revealed the possibility in experimental realization of 
this phenomenon by using silver-germanium multilayered structures, which will stimulate 
experimental efforts to fabricate such anisotropic ENZ metamaterials and characterize this 
exotic optical phenomenon. 

This work was partially supported by the Department of Mechanical and Aerospace 
Engineering, the Materials Research Center, the Intelligent Systems Center, and the Energy 
Research and Development Center at Missouri S&T, the University of Missouri Research 
Board, and the Ralph E. Powe Junior Faculty Enhancement Award. S. Feng was supported 
by NAVAIR's ILIR program. The authors acknowledge Y. He for his useful discussions 
about this work. 



5 



REFERENCES 



X S. Enoch, G. Tayeb, P. Sabouroux, N. Guerin, and P. Vincent, Phys. Rev. Lett. 89, 213902 
(2002). 

2 M. Silveirinha and N. Engheta, Phys. Rev. Lett. 97, 157403 (2006). 
3 M. G. Silveirinha and N. Engheta, Phys. Rev. B 76, 245109 (2007). 
4 R. Liu, Q. Cheng, T. Hand, J. J. Mock, T. J. Cui, S. A. Cummer, and D. R. Smith, Phys. 
Rev. Lett. 100, 023903 (2008). 

5 B. Edwards, A. Alu, M. E. Young, M. Silveirinha, and N. Engheta, Phys. Rev. Lett. 100, 
033903 (2008). 

6 A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, Phys. Rev. B 75, 155410 (2007). 
7 Q. Cheng, W. X. Jiang, and T. J. Cui, Phys. Rev. Lett. 108, 213903 (2012). 
8 A. Alu and N. Engheta, Phys. Rev. E 72, 016623 (2005). 

9 D. Rainwater, A. Kerkhoff, K. Melin, J. C. Soric, G. Moreno, and A. Alu, New J. Phys. 
14, 013054 (2012). 
10 J. B. Pendry, D. Schurig, and D. R. Smith, Science 312, 1780 (2006). 
n D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. 

Smith, Science 314, 977 (2006). 
12 S. Feng, Phys. Rev. Lett. 108, 193904 (2012). 
13 D. J. Bergman and D. Stroud, Solid State Physics 46, 147 (1992). 
14 X. Yang, J. Yao, J. Rho, X. Yin, and X. Zhang, Nat. Photonics 6, 450 (2012). 
15 Y. He, S. He, and X. Yang, Opt. Lett. 37, 2907 (2012). 
16 S. Adachi, Phys. Rev. B 38, 12966 (1988). 

17 R. J. Pollard, A. Murphy, W. R. Hendren, P. R. Evans, R. Atkinson, G. A. Wurtz, A. V. 

Zayats, and V. A. Podolskiy, Phys. Rev. Lett. 102, 127405 (2009). 
18 A. A. Orlov, P. M. Voroshilov, P. A. Belov, and Y. S. Kivshar, Phys. Rev. B 84, 045424 

(2011). 



6 



FIGURE CAPTIONS 



FIG. 1. (Color online) The IFCs of the anisotropic ENZ metamaterial with different ma- 
terial losses, and the simulation results of the loss enhanced transmission in the anisotropic 
ENZ metamaterials. For the anisotropic ENZ metamaterial, the IFC of k x ~ Re(k y ) is plot- 
ted as red solid curve, and the corresponding IFC of k x ~ Im(/c y ) is plotted as red dashed 
curve. The IFC of air is plotted as blue solid curve. The Gaussian beam is TM polarized 
and the distribution of the real part of H z component is plotted. The incoming Gaussian 
beam is denoted by the wave vector k- in and the Poynting vector S- in) with different angles 
of incidence, including 15°, 30°, 45°, and 60°. The propagating Gaussian beam inside the 
anisotropic ENZ metamaterial is denoted by the wave vector fc ENZ and the Poynting vector 

FIG. 2. (Color online) The simulation results of the loss enhanced transmission in 
realistic silver-germanium multilayered stack, placed on a germanium substrate. The shaded 
area represents the multilayer stacking in the x-direction. The incoming Gaussian beam is 
TM polarized and the distribution of the real part of H z component is shown at the ENZ 
frequency 193.4 THz (1.55 fim). The propagation direction of the Gaussian beam is denoted 
by hollowed arrows. Two different angles of incident, 5° and 15°, are considered with respect 
to three different damping rates, including 7, 37, and 97. 
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FIG. 1. Lei Sun, Simin Feng, and Xiaodong Yang 
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FIG. 2. Lei Sun, Simin Feng, and Xiaodong Yang 
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